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AxiSTIl/vCT 


The  forces  on  a  slender  body  of  revolution 
hcverin}-  under  waves  are  deter^inca.  It  is  assumed 
in  the  analysis  that  tna  viavelength  is  of  the  sane 
order  of  maf'.nitude  as  the  body  radius,  i.e.  kPsOd)  , 
where  k  is  the  wave  number  and  P  is  the  bodv 
radius.  Two  cases  are  considereo:  beam  seas  anu  head 
seas.  It  is  shown  that  for  head  seas  conventional 
slander  Dody  theory  is  inadaiuate,  and,  for  this  case, 
a  pcneralization  of  slender  body  theory  is  developed. 
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FORCES  ON  A  HOVERING  SLENDER  SODY  OF  REVOLUTION 
SUBMERGED  UNDER  WAVES  OF  MODERATE  WAVELENGTH 

I.  INTRODUCTION 

Wh«n  a  slander  submerged  body  is  hovering  in  a  velocity 
field  caused  by  a  sinusoidal  wave  train  of  large  wavelength •  the 
forces  on  the  body  can  be  calculated  using  conventional  slender 
body  theory.  Such  a  calculation  has  been  performed,  for  example, 
by  Kaplan  and  Hu  Cl].  In  this  approach  the  orbital  velocity 
potential  associated  with  the  undisturbed  waves  is  expanded  in 
the  neighbo.rhood  of  the  body«  in  which  case  the  problem  reduces 
to  solving  a  two*dinensional  potential  problem  in  the  cross  flow, 
and  the  force  per  unit  length  can  be  determined  in  terms  of  the 
added  mass  of  the  section.  An  implicit  assumption  which  is  made 
in  conventional  slender  body  theory  is  that  the  longitudinal 
variations  are  small  in  coa^rison  with  the  lateral  variations. 

It  is  clear  that  for  a  long  slender  body  in  long  waves  this  as* 
sumption  is  satisfactory  except,  perhaps,  locally  near  a  blunt 
end  or  at  the  beginning  of  an  appendage.  Local  aberrations  can¬ 
not,  of  course,  affect  the  total  force  to  any  great  extent. 

A  different  situation  prevails  when  the  wavelength  is 
of  'ttie  sams  order  of  magnitude  as  'the  dimensions  of  the  cross- 
section  and  ‘ttie  waves  are  not  beam-on.  In  this  case,  longitudi¬ 
nal  variations  become  of  the  same  order  of  magnitude  as  lateral 
variations  and  it  is  necessary,  in  calculating  the  body  potential. 


a. 


to  take  this  into  account*  A  similar  situation  would  prevail  if 
one  were  to  caloulats  the  flow  about  a  slender  corrugated  body 
with  no  waves  if  the  wavelength  of  the  corrugations  were  of  the 
same  order  as  the  cross-sectimial  dimensions* 

In  the  present  report  we  will  deal  with  the  case  of  a 
slender  body  of  revolution t  and  we  will  consistently  make  the 
assumption  that  the  wavelength  is  of  the  same  order  of  magnitude 
as  the  body  radius,  i«e*  kR  «  0(1)  ,  in  which  case,  since  the 

body  is  slender,  kf  1  •  wfe  will  then  seek  the  leading  term 

in  the  expression  for  the  force  consistent  with  tnis  oi*dering 
hypothesis*  In  the  course  of  the  analysis,  body  motions  will  be 
included*  However,  these  motions  are  assumed  to  be  wave-induced, 
and,  as  a  consaquancs,  thair  oscillatozy  part  will  be  at  most  of 
the  order  of  the  wave  orbital  motions;  in  other  words,  an  upper 
bound  on  the  oscillatory  part  of  the  body  motion  is  given  by  the 
motion  of  a  particle*  This  fact  will  aid  us  in  simplifying  the 
expressions  that  are  obtained*  TWo  cases  will  be  considered:  beam 
seas  and  head  seas* 

lbs  velocity  potential  of  the  waves  with  wavelength  X  , 
amplitude  a  ,  and  wave  propa^tion  speed  c  is  given  by: 

where  k  •  2'^^  •  g/c^  and  the .real  part  is  to  be  taken*  The 
coordinates  (x,y,s)  are  fixed  with  respect  to  an  inertial  system, 
with  s  pointing  upward  and  x  dirsotsd  along  the  axis  of  the 
undisturbed  body*  Thm  angle  is  the  angle  between  the  x  axis 
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and  the  wave  create  (see  Figure  I*).  In  tema  of  a  coordinate 
system  fixed  with  respect  to  the  body,  the  wave  potential  becomes 


-  ac6 

where  <  j  •  > 

body  displacements* 
by  (u,v,w). 


6  * 

represent  the  (x,y,x)  components  of  the 
Their  time  rates  of  change  will  be  denoted 


II .  BEAM  SEAS 


For  beam  seas  \  *  T /2  and  Equation  (2)  reduces  to 

^  -  acs'"”  a'"  (9) 


There  is  no  x  variation  in  the  wave  potential  for  beam  seas, 
and  consequently  none  in  the  wave  orbital  velocities*  Furthermore, 
the  X  component  of  orbital  velocity  is  aero*  Because  of  this, 
the  longitudinal  variations  are  indeed  small  compared  with  the 
lateral  variations,  and,  if  the  body  is  slender,  it  is  permissable 
to  use  conventional  slender  body  theory*  The  body  potential, 
therefore,  satisfies  the  two  dimensimal  Laplace  equation  in  the 
cross  flow.  Let  us  transform  to  polar  coordinates  in  the  cross 
flow: 


y  «  r  cos  S 
z  «  r  sin  @ 


(4) 


The  body  potential  can  be  determined  by  using  the  circle  theorem 
[2!]*  The  wave  potential  plus  body  potential  satisfies  the  condi- 
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tion  of  no  noTMl  flow  on  tho  oirelo  v  •  R  •  and  is  givsn  by 

+  coakC^fios^  <5) 

whom  tho  rssl  part  has  boon  takan*  To  this  aost  bs  addad  the 
potantial  induoad  by  tha  body  aotions  w  and  w  (saa  Figura  1.). 
Tha  eoaplata  potantial  is  thus  givsn  by 


aca  cos  k  ^ 

+  (f  vva) 

CO*  © 

r  r 

Tha  radial  and  oireuafarantial  valoeitias  avaluatad  at 
the  body  (r  •  R)  arat 

Cos[«-k'Rcea9 -k«^fkct3 


-  vsin  9  <kwcoso 


V— 


f  CO*  9  f  W  S(n  9 


Tha  tiaa  rata  of  ohanfa  of  thm  vslooity  potential  avaluatad  on 
tha  body  i*  givsn  by 

^  -  Ji  ock «-" 

+  wcosk  CRcoa#4-e|-tf|)3.  ^ 

Hsrs  ths  tarns  proportional  to  v  and  w  havs  bssn  ignorad* 
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Th«se  ignored  terns  give  rise  to  standard  added  nass  terns  in  the 
force  which  can  be  added  in  at  the  and  of  the  analysis* 

The  pressure  on  the  body  is  obtained  fron  the  (two* 
dinensional)  Bernoulli  equation  in  moving  coordinates: 


(10) 


Transfoming  to  polar  coordinates  with  the  aid  of  Equation  (4), 
it  is  found  that 


(11) 


However «  we  will  consistently  ignore  squares  of  the  body  motions t 

Of 

and  since  is  of  the  order  of  the  body  motions  according  to 

Equation  (&)•  the  prsssure  sli^lifies  to 


^  ^  -  w  cos  o  -  i  ^ 

f  r  T©  1  v-r  '5®  f 


(12) 


The  vertical  (heaving)  force  per  unit  length  is  given 


IF 

4»  J  ' 


&*n  $  'R  cLo 


(13) 


Equations  (7*9)  must  be  substituted  into  Equation  (12) 
and»  thence 9  into  Equation  (13)*  .In  so  doingt  all  quadratic  terms 
in  the  boqy  motions  are  to  be  ignored*  In  calculating  the  term 

which  appears  in  Equation  (12) »  it  will  be  neoessary 
to  square  Equation  (7)^  and  in  so  doing  the  term  cos^(9 ‘‘k1lc«am*k*|^'4kct') 
will  appear*  The  double  angle  formula 
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(14) 


must  b«  used,  and  the  higher  harmonic  term  will  be  ignored*  The 

heaving  force  per  unit  length  is  then  given  by 

f  k^Slr»©  /  ,^\ 

|(v-C.^)e  Sin  K  )  Sir.  (9  0^ 

•  w  (‘Rco&<P  si«i  tp  cl*®] 

O 

0 

"V^iiii^S  €  co6(0 -kUco^®  -k'>^  J 

r  (  .-1<C^-  f  Xk'R.fci.iC? 

iack^Tle^^''  "^Mo^ke  ]e  s.;n  ®  ^  ss) 


k  • 

+  2w  Jcos®  sin©  Cps(©  -klicos©  •  k<v^'|.k6'T)  d® 
r .  1  kiisi«®  ,  ^  .  \  '  "1 

"  2yj  sin  ©  e,  Ccs^®  -kKCi?S»  -k'vj  )  01® 

9  “* 


(15) 


There  are  five  different  integrals  in  Equation  (15)  an< 
each  of  them  is  determined  in  the  Appendix*  The  final  result  is 


where  is.  a  modified  Dessel  function  of  the  first  kind* 

Tlie  first  term  in  this  equation  is  a  suction  force*  It 
is  of  second  order  in  wave  amplitude*  and  attenuates  twice  as 
rapidly  with  depth  h  as  the  second  term  which  is  oscillatory* 
Nevertheless*  for  certain  operations  it  is  important  to  retain 


the  suction  force  because  it  pez'sistently  acts  upward  and  tends 

to  make  the  body  rise.  Motice  that  explicit  dependence  on  the 

1 

body  velocities  (v,w)  has  disappeared.  The  result  still  de¬ 
pends  on  body  motions,  however,  through  the  displacements  )• 

A.  An  Alternative  Derivation 

The  derivation  of  the  heaving  force  given  above  is  per¬ 
fectly  straightforward,  we  will  now  present  an  alternative  deri¬ 
vation  from  a  different  point  of  view  which  is  actually  a  prologue 
to  the  derivation  of  the  heaving  force  for  a  body  in  a  random  sea. 

Consider  the  body  displacements  j  ^  and  the  body 
velocities  v,w  .  From  the  equations  of  motion  of  the  body  we 
can  expect  that  the  body  motiwis  will  consist  of  a  drifting  motion 
which  varies  slowly  with  time  and  which  is  due,  primarily,  to  the 
suction  force,  plus  an  oscillatory  motion  which  varies  rapidly 
with  time  and  which  is  due,  primarily,  to  the  oscillatory  force. 
Hence,  let 

I 


■/  »  /c  '  ^  ' 


(17) 


W«  Wc’t'  w' 

where  the  quantities  with  subscript  zero  denote  the  drifting  motims, 
and  the  primed  quantities  denote  the  rapidly  oscillating  motions. 

The  velocity  potential  and  pressure  can  be  similarly  separated, 
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If  the  relations  between  all  quantities  were  linear  then  we  could 
state  unequivocally  that  the  oscillating  input  caused  the  oscillat¬ 
ing  output  euid  the  drifting  input  caused  the  drifting  output. 
However,  Bernoulli's  equation  for  the  pressure  is  quadratic  with 
respect  to  the  velocity  potential,  cund  even  the  potential  itself. 
Equation  (6),  involves  the  body  motions  )  in  a  nonlinear 

way.  It  is  clear,  therefore,  that  the  drifting  part  of  the  input 
can  affect  the  oscillatory  part  of  the  output  and  vice  versa. 


Substituting  Equations  (17,18)  into  Equation  (6)  there  is  obtained 

I  -U  (ir.  -  rjiirt  9)  . 

^  c.<?5k  f  C660  +  ft*  “CT) 


•f-CVc+  y'^  C0&5?  W') $in© 


(19) 


Similarly  substituting  into  Bernoulli's  equation.  Equation  (12), 
there  is  obtained 


4e  will  retain  terms  which  may  be  of  the  order  of  the  wave  ampli¬ 
tude  squared  and  ignore  higher  order  terms.  In  so  doing,  we  will 
assume  that  the  oscillatory  motions  are,  at  most,  of  the  order  of 
the  wave  amplitude  .  This  is  a  reasonable  assumption  in  view  of 
the  fact  that  the  motions  are  solely  wave-induced.  Thus  in  Equa¬ 
tion  (19),  we  expand  the  exponential  and  trigonometric  terms 
treating  and  ai'  as  small  quantities,  and  retain  only  linear 
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terms  in  the  expanalcm. 


I  (  Lf  ^ 

^  4.  ^  »  ccc.  \  e  c,o6  « {.  r tofe  O  +  -ct) 

4-k^'&  W  ^1*  CoS®  T-'^o-C.'t) 

Je  **^''’®^STnk  Lrtos® 

'•”  "“  s**'©)  ,  /  x\ 

^  cos  Ic  (.“?■  ®  ■*■ 

f  ^‘”®V«>sic[^-cps0  VV|. - 

+  C'^0+  v‘)  ^  cos  S  -sfw,  t  W*')  ~  sin  ® 


-  k  '^1  e 
f  e. 


(21) 


Now  average  Equation  (21)  over  one  period  and  denote  the  averag¬ 
ing  with  a  bar.  Clearly*  the  average  of  any  primed  quantity  is 
zero.  Drifting  quantities  will  be  assumed  to  be  almost  constant 
over  one  period  and,  therefore,  unaffected  by  the  averaging  pro¬ 
cess,  They  are  said  to  sift  through  the  averaging  integral, 

-kCk-^,-»'i‘ne)^ -  - - - 

acKe  ^i'cesk  Crcos  ©  ^  ^  ®-  + 

'•■Q.ckC  T*'‘‘®^[^'c«sk(-^coSS+^j,-c.'t'j-  y  sin  kL^CoS® 

+  Vo  B.  cosd  rw,  ^  Sin®  (22) 

r  r 

Subtracting  Equatim  (22)  from  Equation  (21)  we  obtain  the  follow¬ 
ing  equation  for  : 
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acle'^^*^  ^,-rSmo  j, 

+  k£',*'’'^‘'’'^’'"‘'''»Uk  Crcoso  +  <^.-«) 

fct  .kCW-Jf.-rtln©^  ',■  T 

-  lO  ft  £osi<ir<:.o&©f''|»- fi.t ) 

,  •W(U*^o*rWnS) 

”  K ®  sin  k  C<*c-o^.C5  4'y  , -ct ) 

,  -k  -^j-rsln®)  ~ 

+  k  o  ft-  S^o  k  Cr  c.o<>  <9  f  ♦!,  -  ct) 

-UCh- i,-  S! 

»  ^  Co6  W  {  '~COS  &  4'<‘ 

^ ^  ^S3Sk  C~ir  t-cs  g 

,  / '  -k  C*^  ^  si t-»e)  .  /  9 '•  ^  ,  77 

•“  K^A  r  cosk  V.  'i» '“^'^ ) 

I  -kCh-X*  ^  ,  I  /ic^  ,\ 

-  k^e  Sm  k  (J^cob®  f  *-J,-ct) 

^  ^*1  ft-  Oink  (_-p  cOb  ® 

+  \/  ^  Coi  &  t  Vi/  ^  Sin  O 

r  r 


(23) 


Th«  first  and  sixth  tsrs  in  the  bracket  must  be  retained 
but  all  other  terms  can  be  ignored  since  they  give  rise  only  to 
higher  harmonics.  Whence* 

I  r 

J«ac-|d  CoiK\^rCo»6>rV»“ft'^)  + 


+  e  7’^"^®VosW(^C3s©-'y..c.t)^ 

+•  V  -p  cos  ®  ’♦'  vJ  ^  ® 

r  r 


(24) 


To  detsxsdne  the  mean  pressure*  average  Equati<m  (20) 


over  one  period. 


•k  (  V  Sin  d  -  wV  rs  0  )  -p 


(25) 
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The  term  can  be  obtained  from  Equatlm  (22).  In  perform¬ 

ing  the  differentiaticmt  it  will  be  assumed  that  i. 
filter  out  of  the  averaging.  However,  since  the  first  four  terms 
are  already  of  second  order  in  the  wave  amplitude,  and  since 
and  are  wave  induced,  any  teimis  which  arise  from  the  filter- 

Ing  of  and  will  be  of  higher  order  and  can  be  ignored. 

Terus  which  arise  from  the  filtering  of  c  ,  on  the  other  hand:  , 
must  be  retained. 

The  terms  proportional  to  v^  and  w^  have  bean  ignored,  consist¬ 
ent  with  our  policy  of  not  considering  added  mass  terms. 

Ignoring  terms  of  the  order  of  the  body  motion  squared, 


the  mean  pressure  beoomes 


J-T.  -  (/»ine  -w'coi  ±  ^  ^ 


(27) 


To  obtain  p  subtract  Equation  (27)  from  Equation  (20), 
If  all  terms  of  seoond  order  in  the  body  motion  and  all  highar 
harmonics  are  ignored,  there  is  obtained 

j'P'-  (v.  5m  9 

When  substituting  into  Equation  (28)  it  is  pexsdss- 
ible  to  ignore  the  first  four  terms  in  Equation  (22)  because  they 
give  rise  to  terms  of  higher  order  in  the  pressure.  Then,  upon 
evaluating  at  the  body,  Equatlm  (28)  simplifies  to 
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Sirt#  »  W,Co&  s)~ 


(29) 


Upon  substituting  Equations  (24,26)  into  Equation  (27) »  lotting 
r  ■  R  and  intsgrating  with  respsct  to  according  to  Equation 
(13),  ths  Man  vortical  forco  par  unit  length  bocoMS  aftor  high¬ 


er  ordor  torns  aro  ignorodt 

A -7  -  .  .  t.  ^ 


(30) 


The  fluctuating  vortical  force  per  unit  length  is  obtained  by  in¬ 
tegrating  Equation  (29)  with  respect  to  ,  and  there  is  obtained: 

ii-  -  i  e  ^  e  **  k  (ct  -  ^  ^  3 j 

4x 

It  can  easily  be  seen  that  Equation  (31)  is  identical  to  the 
second  tern  of  Equation  (19)  with  and  replacing  f*!  and 
^  respectively*  The  first  tern  of  Equation  (30)  is  the  sane 
as  the  first  tern  of  Equation  (16)  with  J,,  replacing  ^  •  The 

second  and  third  terns  of  Equation  (30)  are  inplioit  in  Equation 
(16)  also,  for  if  we  let  -f  and  ^  /  ^  and 

then  expand  Equation  (16)  retaining  terns  of  the  order  of  the 
square  of  the  wave  anplituds,  the  seoond  and  third  terns  of  Equa^ 
ti«i  (30)  are  rsproduoed  idsntieally  once  higher  harnonics  are 
ignored.  The  oscillatory  force  given  by  Equation  (31)  is  identi¬ 
cal  with  the  result  which  would  have  been  stained  if  the  Imig 
wave  approximation  were  nade  at  the  outset  of  the  analysis  (see 
CD).  Zt  can  be  shown  that  all  other  first  order  forces  and 
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Boaents  are,  likewise,  identical  to  their  long  wave  limits  in 

beam  seas* 

If  the  eseillatory  part  of  the  BOtion  is  such  that  the 
body  Boves  like  a  water  particle  (which  is  very  nearly  the  case). 


then 


I  -w  1/4.  \ 

•v)  «  a  e  CoS  K  Cct  -  ) 


(32) 


Substituting  into  Equation  (30),  there  is  obtained  the  following 
expression  for  thm  aean  suction  force  per  unit  length ■ 

- Ij  « 

th.  univara.l  funetiGn  —  I  J  1. 


(33) 


tabulated  in  Table  X* 


II«A*  RANDOM  SEAS 


Suppose  the  wave  potential,  as  given  by  Equation  (1) 
for  a  sinusoidal  wave  train,  is  replaced  by  a  function  of  a  ran- 
doB  variable  with  spectruB  A*’(fa9)  •  Using  the  formulation  of 

Pierson  CH],  the  coi^lete  potential.  Equation  (6),  is  replaced  by 


- 1*)} f  jKcoi.®  4. »In « 

Now  assuBe  that  all  quantities  can  be  written  in  ^e  form 


-kCli-  S-Sin^) 


(31) 


(  )  •  (  )o  ♦  (  ) 


05) 


*This  result  has  been  <Atained  independently  by  Ogilvie  C3]* 
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where  the  subscript  zero  represents  the  mean  value »  and  is  a  slow¬ 
ly  varying  function  of  time*  Quantities  with  a  prime  are  fluctuat¬ 
ing  random  quauitltles  with  zero  mean*  The  similarity  between  this 
notation  and  the  notation  used  In  the  preceding  section  will  be 
noted* 

Continuing  the  parallel  with  the  preceding  section,  we 

will  operate  on  the  potential  and  the  pressure  with  the  average: 

T 

J —  \  (36) 

This  average  replaces  the  average  over  one  period  which 
was  used  In  the  preceding  section.  By  the  limit  T-^^  we  meam 
that  T  Is  large  enough  so  that  a  meaningful  sample  of  any  random 
quantity  Is  included,  but  small  enough  so  that  mean  quantities  are 
essentially  constant  over  this  time  Interval*  Thus,  quantities 
with  subscript  zero  will  filter  through  the  operation  represented 
by  Equation  (36) ,  just  as  they  did  for  the  average  over  one  period 
in  the  pezdodic  case.  The  analysis  Is  Identical  In  every  respect 
to  the  analysis  performed  for  periodic  fluctuations,  provided  only 
that  we  interpret  the  bar  to  be  the  etatlstlcal  average  represent¬ 
ed  by  EquatiOT  (36). 

Suppose  the  fluctuating  part  of  the  notion  to  be  such 
that  the  body  follows  a  water  particle.  Then  the  equation  corre¬ 
sponding  to  Equation  (33)  is 


(37) 


15. 


llie  fluctuating  force  per  unit  lengti)  ie 

$  J'lc Ve”** * Vm  l<  -  '^o)  Am 


(38) 


If  the  wave  spectrun  is  a  Neumann  spectrum  for  a  fully  developed 


sea  (Reference  C4j) 


Then 


iZi  .ZillQ 

in  0  J« 


A  C«^)  -  I  e.  (39) 

Then,  letting  c.*  ,  k«  ,  Equation  (37)  becomes 

If  the  waves  are  long  with  respect  to  the  body  radius •  then 

1+  '  t  and  the  a>  integral  can  be  evaluated 

«''actly«  The  integrated  force  can  then  be  evaluated  to  be 


(41) 


where  is  the  voluastric  polar  monent  of  inertia.  If  the  waves 
are  not  long^  but  the  body  is  deeply  submerged  ** 
the  4>  integral  can  be  evaluated  approaciaately  by  the  method  of 
steepest  descents  CS]*  The  result  is 

f  •*  (M) 


where  the  integral  is  taken  over  the  length  of  the  body. 
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III.  HEAD  SEAS 

For  head  seas  *  0  •  and  Equation  (2)  reduces  to 

ace  e  (43) 

There  is  an  explicit  x  variation  in  the  wave  potential  for  head 
seas*  but  no  y  variation.  C^sequently,  there  are  x  varia¬ 
tions  in  the  orbital  velocities*  but  the  y-coaponent  of  ozbital 
velocity  is  zero.  The  wavelengtii  of  the  orbital  velocity  will 
be  assuxed  to  be  of  the  sane  order  as  a  cross-sectional  diaension* 
i.e.  )cK  «  0(1)  *  and*  as  a  consequence*  conventional  slander 

body  theory  is  not  applicable*  For  convenience*  let 

21  *  X'  J  (44J 


and  transfora  to  polar  coordinates  according  to  Equation  (I). 
The  wave  potential  then  bacmMSs 


(45) 


Expanding  the  wave  potential  as  a  Fourier  series  in  d  *  there 
is  obtained 

where  is  a  aodified  Bessel  function  of  the  first  kind. 

The  body  potential  nust  satisfy  the  three  diaensicaial 
Laplace  equation.  In  cylindrical  coordinates*  tills  potential  can 
be  expressed  in  terns  of  line  sources  and  higgler  order  line  singu- 
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lariti«s«  as  follows: 

>  I  m.  I 


J  [Cx-.y.  rT ^  "  4-  ^  ^  i  r^.-x.)V  r^]" 


(47) 


The  potential  will  be  used  to  cancel  the  normal  flow  at  the 

body*  r  •  R(x)  •  It  is*  likewise*  necessary  to  add  another 
potential*  •  to  account  for  the  body  motions.  However*  the 

body  motions  induce  velocities  which  are  slowly  varying  in  the  x 
direction*  and*  consequently*  an  ordinary  slender  body  representa¬ 
tion  will  be  adequate  for  this  potential.  In  other  words*  the 
potential  id  also  represented  by  an  equation  of  the  form 
of  (47) «  but  we  nay  take  the  limit  rrb>0  without  regard  to  the 
rapidity  of  the  variation  of  f  with  respect  to  x,  ,  This 
limit  yields  (see  Reference  fS]): 

f  a  -  —  +  qCx)+I(f2."^ - - -  (‘»8) 


(48) 


where 


(49) 


The  expression  for  can  also  be  simplified  by  taking  the  limit 
r-»  0  *  but*  because  this  potential  must  oanoal  the  wave  orbital 

velocities*  the  strength  of  tits  singularities  will  vary  rapidly 
with  X|  •  Uenoe*  let 

(SO) 


(SO) 


where  F  is  a  slowly  varying  f\motion  of  x  ,  ,  Equation  (47) 
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bets 


^  '  “  LC'-X')  f-r  ] 


ikx, 

dki 

o^rr^ 


(51) 


We  nust  now  ta)ce  the  limit  ,  which  expresses  the  fact  that 

the  wavelength  is  small  compared  with  the  body  lengthy  i.e*  ki  >  1. 
According  to  a  lemma  of  Riemonn-Lebesque •  Reference  [y3»  integrals 
of  the  form  found  in  Equation  (SI)  will  vanish  as  k’^eo  ,  pro¬ 
vided  the  non-oscillatory  port  of  the  integrand  has  bounded  varia¬ 
tion.  However^  if  we  simultaneously  take  the  limit  r  0  ^  it 

is  clear  that  all  the  integrands  will  become  singular  in  tiie 
neighborhood  of  x,  ex  •  Hence »  only  that  port  of  tiie  integrand 
in  the  immediate  vicinity  of  the  singularity  will  contribute  to 
the  integral.  We  may.  therefore,  evaluate  ^ 

and  bring  Fu  outside  the  integral  sign.  Furthermore,  we  may 
let  the  limits  of  integration  extend  from  x,*-o^  to  Xia-t-od 
sinos  tise  results  of  so  doing  will  contribute  a  negligible  amount 
to  the  value  of  ^e  integral  by  virtue  of  the  Riemonn-Lsbesque 

lemma.  In  aocordonoe  with  these  concepts.  Equation  (51)  beoomss 

'  «»  lluf. 

if  . 

(52) 


The  intsgrals  oon  be  expressed  in^  terms  of  modified  Bessel  funo- 
tions  of  the  second  kind,  with  the  result  that 

V  rco+x) 


19. 


In  the  neighborhood  of  the  body*  the  quantity  (kr)  which  ap¬ 
pears  as  the  argumnt  of  the  Bessel  functions*  is  postulated  to 
be  of  order  unity*  and*  as  a  consequenos*  Equation  (S3)  cannot 
be  simplified  any  further.  The  complete  potential  is  given  as 
the  sum 


(Sl>> 


where  expressions  for  •  'fm  given  by  Equations 

(46*48*54)  respectively. 

The  boundary  condition  for  a  three-dimensional  body  of 
revolution*  including  the  body  motions  surge*  heave*  and  pitch* 
has  been  derived  by  Cuthbert  and  Kerr  in  Reference  CfJ: 

;fx]  ^  *(w-  ^)t)cost«9-\)+  Jr;  r-X  <55) 

where  the  subscripts  x*r  denote  partial  derivatives.  By  eem» 
paring  the  second  term  on  the  left  with  the  first  term  cn  the 
light*  it  is  clear  that  the  second  term  on  the  left  can  be  neglect¬ 
ed  for  a  slender  body*  and  the  bcnindary  condition  then  sii^lifies 
to 

Differentiating  Equations  (46*48*S4)  with  respect  to  x  and  r  * 
and  substituting  into  Equation  (56)*  there  is  obtained 
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,  1  K«)c  (w'"'*‘ff  it  (*“'f-“)(i  W-tiK-fr'Vf (.'Ol  r 

\  Ofy_C^  PMcosh  (g-T^) 

’  TO  4"  r*/.- .  I  \ 


(57) 


^  rcH.o 

+ 1  l<K  Cy)d  K>  Ck^) + ik iT  zl  ■  \  \ 

Con8id€r  the  t«nii8  corresponding  to  n  «  0  » 

[u.  V  ocilce ^  (Ke‘'‘ ')]  ^ 

,  -i.(y.-J)ftkZ  i  .  -  LV k’  fl.tt^  <s«) 

=  \<tuce  IfllkR)  i 


'  tU**-  CLC,t> 

~  w 

+  kl  (kR) 

O'*  •  ^  i  "P^ 

We  equate  the  non»oscillatory  parts* 

To  the  first  order  in  hody  radius  we  obtain 

f.Cx>  .  - 


(59) 


u 


(50) 


We  now  equate  the  oscillatory  parts* 

ixL 

i;(ku')+ ikF.e*'''  k.'  (kR^ 


(51) 
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If  XR  ■  0(1)  •  th«nt  to  tho  firet  ordor  in  body  rndiu*  tra  db- 


t4d.n 


Ikx 


f  (62) 

Now  consider  tho  non-oscillatory  terms  corresponding  to  n  «  1  . 

To  the  first  order  in  bo<!^  radius  we  obtain 


f,-  n'w 


(63) 


Finally^  consider  the  oscillatory  terns  corresponding  to  n  >  1  . 

To  the  first  order  in  body  radius  we  obtain 

,'k«  F.(,)  -  - 

Substituting  back  into  the  potential*  we  obtain  after  sone  com¬ 
bining  and  siaplifying 

J  .  K’.(W) 

Where  g(x)  «  -2u <l)r,  •  “** 

C  s I 

^  ti 

*1  Y\}0 

The  partial  derivatives  of  (j>  with  respect  to  r»  dP  *  » 
and  t  evaluated  at  the  body  r  »  R  *  are  given  respectively  by 


(65) 


^  -U-^-Wcos(^-^) 

Or  ix  ^ 


(66) 
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7^  WV*  4- 

.  -k(w-i:)*ikS  «  ^  ' 

K'»  cui') 


(67) 


0^  ot&ta 
Ox 


K*.  CklO 

Z.  If'.futo  *» 


¥- 


+  Tc-|r(tf-)^sCs-50n'^^ 

cosnC^-!^^ 

k;  ckTi'i 


(69) 


wh«r«  th«  Wponskian  identity « 


I 


(70) 


hM  ba«n  aaployad  it««  Rafamnca  C93).  In  additions  highar  ordar 
tama  in  bady  radiua  hava  baan  dvoppad  in  the  axpraaaien  for  -^2- , 

oiC 

and  tanns  which  contributa  to  tha  addad  aasa  forea  hav*  baan  ig* 
novad  in  tha  axpraaaion  for  • 

Tha  praaaura  on  tha  body  ia  obtainad  froai  lha  (thraa* 
diaanaional)  Bamoalli  aquation  in  aoving  ooordinataa.  E^prasaad 
in  taraa  of  oylindrioal  ooordinataa  this  ia 

V'  T-  ^  ®) 

-(«-+ ^  -  X  [(t-  -*•  (^y] 


(71) 


But  ia  of  tha  ordar  of  tha  body  aotiona  aeeerding  to  Eqtui- 

tion  (66) •  and  wa  will  oonais'tently  ignora  squaraa  of  tha  body 
aotiona*  Uanoa  tha  praaaura  si^lifiaa  to 


(72) 


^  “(w -<V») -p  ^COi  9 -(a  V  ratio  ©)  ^ 

-Xl^r 

Tha  vortical  foroa  par  unit  langtii  ia  givan  by  Equa^on  (13)* 
Aftar  taking  tha  raal  parts  Equationa  (67*66)  naat  ba  aubatitutad 
into  Equation  (72)  and  thanoa  into  Equation  (13)*  Zn  ao  doings 
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V 


«11  qtt«dr«tie  tcras  in  body  aotions  ar«  to  bo  Ignorod,  and 
ttim  ofthogenal  pveportioa  of  ^o  trigcnoaotrio  tonia  in  B  aay 
bo  Mod*  Zft  at  tho  oaM  tiao,  wo  ignoro  highor  ordor  tbrao  in 
body  radiM^  wo  aro  finally  lad  to  tho  following  oxproaoien  for 
tho  baaving  foroo  par  unit  long-tiit 

jZ  f-  I  r 

fix  n  4- <;(kttK’,»i(ifK)L  C'®)*  J 

Ki(W  “  <■ 


Zn  arriving  at  Equation  (73)  it  has  boon  asauaa.d  that  tho  aurgo 
valeoity  u  ia  of  tho  aano  ordor  of  Mgnitudo  aa  tho  haaviag 
volocdty  «  •  in  which  eaaa  tho  spooiflc  dopondoneo  of  tiio  vorti- 

oal  foroo  on  u  ia  nogligiblo  for  a  alondor  bo4y* 

Equation  (73)  can  bo  ainplifiod  ovon  furthor  by  talcing 
into  aooount  tttm  ordor  of  nagnitudo  of  tho  bo4y  aotiona  (whiidk 
ara  wava  induoad)  i  and  alao  by  fooMing  attantion  on  tho  total 
foroa  inataad  of  tbm  foroo  par  unit  langth* 

Conaidar  tho  body  aotiona «  and  lot  oadk  notion  ba  rapro- 
aantod  by  a  alowly  varying  part  plM  a  rapidly  oaeillating  part* 
in  tho  aaaa  nannor  aa  for  tho  boon  aaa  eaao*  For  hoad  aaaa»  tho 
oaeillating  part  of  ttM  notion  ^11  oortainly  ba  aaallar  than  tho 
notion  of  a  watar  partiola  and  ao  tha  oacillatory  diaplaoananta 
given  by  Equation  (32)  ropraaont  upper  bounda*  An  ivpar  bonad 
on  tha  oaeillatory  valocitioa  will  than  bo 

J 


“Oucke  (ci- 


(7a) 
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Th«  drifting  vnlooitics  which  arc  wava-induead  ara 

clearly  much  aaallar  than  the  wave  apaad  e  t  hanoa  the 
drifting  tarwa  which  appear  in  ‘^a  coa  kj(.  tana  of  Equaticn 
(73)  ara  nagligibla  coaiparad  with  tha  antira  ain  tarw. 

Upon  aatiauiting  tha  oacillatory  valocitiaa  with  tha  aid  of 
Equation  (74)  •  wa  obtain  an  uppar  bound  for  tha  coa  tarn 

of  Equation  (73) 

K; 

which  ia  aaan  to  ba  of  tha  aana  ordmr  of  nagnituda  aa  tha  firat 
tana  of  Equaticn  (73)«  i.a.  0(^  a'‘c^k)  •  On  tha  baaia  of 

this  arguaant»  we  nuat  retain  'Oic  coa  k  JL.  tarn.  Howawar,  the 


(75) 


total  force  ia  tha  integral  over  tha  length  of  tha  body.  Inte¬ 
grating  tha  firat  tarn  of  Equation  (73)  over  tha  length  of  tha 
bodyt  it  ia  aaan  that  thia  tarn  ia  0(^a^e^k  i  )  t  wham  tha 
ordaring  hypothaaia  kR  •  0(1)  has  bean  utilisad.  On  tha  othar 
handy  integrating  Equation  (75)  over  tha  length  of  the  body. 


them  ia  obtained 


A  -v  cotk  (k ♦  ^ 


(76) 


According  to  tha  Rianann-Labaaqua  lanawi  (Rafarenoa  C7])a  as 
k-.oo  tha  intagml  ia  0(l/k)  .  Hanoa*  thm  coa  kX  tarn  is* 

In  affaoty  ^  a*  c^  )  .  Tha  ratio  of  thia  tarn  to  tha  fimt 

term  of  Equation  (73)  is  0(lA£  )  which  nay  ba  naglactad.  For 
a  closed  body  1/k/ (kR)  vanishes  at  tha  and  points  of  tha  inta- 
gralftandy  in  this  cassy  if  tha  body  mdius  has  a  continuous  fimt 


derivative y  it  can  ba  shown y  by  an  integration  by  parts*  that  tha 
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ratio  of  the  coa  kj(  term  to  the  first  term  of  Equation  (73) 
actually  vanishes  as  (l/k£  )^  •  A  similar  argument  can  be 
used  to  show  that  the  oscillatory  parts  of  and  J  may  be 
neglected  in  the  exponentials  of  the  remaining  terms.  Hence » 
the  total  heaving  force  becomes 

Z - 4- ie [s(.rW^(kIl) sin  k 

(77) 


(78) 


CU) 


I 

k;  (kB.) 


(79) 


and  the  integrations  extend  over  the  length  of  the  body.  The 
universal  functions  Ck,  and  are  tabulated  in  Table  I. 

The  result  for  head,  seas  may  be  contrasted  with  that 
for  beam  seas.  In  the  beam  sea  case  the  oscillatory  body  notions 
contribute  terms  of  the  sans  order  of  magnitude  as  the  suction 
force  on  a  stationax*y  bo<l^»  and»  furthermore,  the  body  motions 
tend  to  reduce  the  suction  fores  by  reducing  the  relative  motion 
of  body  to  water.  In  the  head  sea  ease ,  the  body  motions  ars  not 
necessarily  in  phase  with  the  water  motion,  and  hence  do  not 
necessarily  rsducs  the  suction  force.  But,  in  any  event,  the 
effect  of  the  bo<fy  noti<»is  is  of  higher  order,  and  may  be  neglect¬ 
ed.  Because  of  the  relief  afforded  the  suction  force  in  the  beam 
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••a  CM«,  w«  may  axpect  tha  suction  fores  to  bs  snallsr  in  bsam 

ssas  than  in  hsad  ssas. 

In  tha  long  wava  limit  tha  ordaring  hypothasis  is 
KR  ^  1  •  k£,  •  0(1)  .  Tha  analysis  has  baan  carriad  out  by 
Cuthbart  and  Karr  [8]  for  long  wavas,  and  thay  damonstrata  that 
tha  affact  of  body  motions  is  prasant  in  tha  foroa  par  unit 
langth  for  haad  saas  as  it  is  in  lha  prasant  analysis.  Howavar, 
if  in  tha  long  wava  analysis  wa  taka  tha  limit  k-*oo  and 
assuma  that  body  motions  ara  boundad  in  aceoxHlanoa  with  Equation 
(7h),  tha  affact  of  body  motions  on  tits  intagratad  fores  in  haad 
saas  is  nagligibla.  If  in  tha  prasant  analysis  %»a  lat  KR~*>0 
than  1*S  t  Equation  (77)  raduoas  to  tha  long 

wava  solution  for  k  oo  . 

Tha  pitching  mosmnt  can  aasily  ba  obtainad  by  introduc¬ 
ing  an  X  into  tha  intagrand  in  Equation  (77).  Tha  lataral  foroa 
and  yawing  moswnt  srs,  of  eoursa,  ssro  for  haad  saas. 

In  a  similar  fashion •  and  naglseting  tarns  of  tha  ordar 
of  tiia  square  of  ^a  wavs  anplituds*  wa  can  dsriva  tha  following 
sxprsssion  for  tha  surgs  fores 


whsrs 


(80) 


(81) 


Tha  univarsal  function  ci^(lcR)  is  tabulated  in  Table  Z.  In  tha 
long  wavs  limit  — >  1  •  and  Equation  (80)  rsduoss  to  tha 
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eoTMOt  long  wovo  solution  (ooo  Koforonoo  Cll). 

IZI.A.  RANDOM  SEAS 

Argunsnts  sinilar  to  thot  usod  for  boon  seos  lood  to 
tho  follouing  oxprsssions  for  tho  stoody  suction  foroo  in  hood 
0008  whon  tho  wovos  ors  dsseribod  by  o  spoctrun  A^'CCO) 

2.-  ‘*** 

^  • 

Tho  fluotuoting  foroo  is 

2  -  25 1 Va‘M 

(»3) 


Tho  fluotuoting  surgo  foroo  is 

x'-  -^{JFUsiTZ  [s6.w,(iiA)  c^sk(x-^,-c.i)Ax 

(M) 


Using  o  Nouaonn  spoctrun  (Equotion  (39)>»  tho  stoody  suction  foroo 

bOOCMO 

dlctf 


(IS) 


2.*  ^  A-  It?  •*  r{- 

tar  long  wovoo  l.i  ond  tho  *0  intsgrol  oon  bo  ouoluotod  sn¬ 
oot  ly 


2  ^ 

®  4  V  /  tj 


(•«) 
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For  th«  dooply  subaargod  e«8«»  the  method  of  steepest  descents 
yields 


(87) 


TABLE  I. 


T«bl«s  of 

th«  ftmotions 

c^M)  - 

I,  (1WR> 

UR 

1 

- 

.  1  ^  ■ 

uv _ _ 

- 

1 

icR'  k;  i  w 

2.  r  d.. 

kR"  J 

• 

kR 

dl«(kR) 

C(.(kR) 

C^^(kR) 

<3l^(kR) 

0,0 

0.00000 

1.5000 

-1.00000 

1.000 

0.1 

0.00501 

1.5035 

-0.99066 

1.009 

0,2 

0.02013 

1.5211 

-0.97532 

1.026 

0.3 

0.0H563 

1.5569 

-0.96126 

1.050 

O.H 

0.00216 

1.6127 

-0.95051 

1.060 

0,5 

0.13032 

1.6699 

-0.96600 

1.115 

o.e 

0.19113 

1.7900 

-0.96197 

1.156 

0.7 

0.26585 

1.9151 

-0.96662 

1.196 

0.1 

0.35601 

2.0676 

-0.95125 

1.267 

0.9 

0.H6352 

2.2506 

-0.96236 

1.300 

1.0 

0.59066 

2,6660 

-0.97756 

1.356 

1.1 

0.76009 

• 

2.7262 

-0.99689 

1.622 

1.2 

0.91510 

3.0256 

-1.02026 

1.690 
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APPENDIX 


In  Equation  (15)  a  number  of  definite  integrals  must  be 
evaluated.  We  will  show  how  these  integrals  can  be  obtained. 


Consider  the  integral 
CP  sir*  ©•  ^  (a-0 


(Al) 


If  the  substitution  (9  *  Ck  *  f 2  is  fflade»  the  integral  can 

be  reduced  to 

T 

r  ik'K.feSdi. 

2-  VcOSCkfi.  ilcIL 

0 

Which •  according  to  formula  (H)  on  p.  181  of  Reference  [63 •  is 


iT  r, 

Consider  the  integral 


(A2) 


ktlsino  r*  1 I  1  XT  I  (A 

'P  e.  -KT.coi9-K*>^+KC.tJ/aln^  60*®  dl© 

If  the  substitution  &  •  ci  t  ^  /2  is  made*  the  integral  can 


(A3) 


be  reduced  to 


Co5l<(ci‘''»j)  d  sin  IcK 

By  using  the  product  into  sum  formula  for  trigonometric  functions 
this  can  be  expanded  into 

»  ®  “  A 


A2. 

In  ordar  to  ovaluoto  thoto  two  iiitogpals»  oonsldor  tho  following 
idontity  (ooo  Riforonoo  [!])•  ^ 

^  » 
lot  y>»x  y  and  wa  obtain 

-  t/  »^o 

•  -  0  y  < «  (AX) 

By  applying  Equation  (AH),  it  ia  aaan,  finally,  that  tha  intagral 
prasantad  in  Equation  (A3)  xo<hioa8  to 

^IcUcptk 

Similarly,  tha  intagral 

IT 

^ CoS [©-kl^CoS®  -lc4j+kot^  d©  (A6) 

oan  ha  shown  to  ha 

(A7) 


(At) 


can,  similarly,  ba  shown  to  ba 

rkn  [(«r-c)  Smk  (ci"  •  +  W  cosW  (ci  -  '»|) J 


- ^  VV.  sink 

Also  tha  tsrms 

XF  * 

-  w  J  e'^'‘*‘"*Cosk  (11»CP4  •  -  at)  sine  de 


(At) 


